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Theoretical Studies of Solar Lasers and Converters 

A second computer program has been developed for the simulation of an n — C 3 F 7 I iodine laser. This 
computer program is given in the Appendix A and typical output from the computer program is illustrated 
in the Appendix B. 

Chemical Kinetics 

The computer program simulates the chemical kinetics occurring during the operation of an n — C 3 F 7 I 
iodine laser. Letting R denote the radical n — C 3 F 7 , these reactions are summarized as follows: 
Photodissociation of RI and I 2 : 

0.27 fxm 

Rl+hu 1 , R + r 

51 

0.5pm 

I 2 + hi/2 / + I* 

52 

Laser action (stimulated emission) 

r + hu 3 ^r 1 + 2/11/3 

Quenching of I* (reaction rates have units [err? /see]) 


r + r/-9X 

1 

+ 

RI 

r + i 2 ^ 

1 

+ 

h 

i* + r 

1 

+ 

R 

r + r 2 -^» 

1 

+ 

Ra 

r + A 

1 

+ 

I 


RI generation and two body recombinations (reaction rates have unitB [cm 3 / sec]) 


R + r -ii* 

RI 

R + /-^+ 

RI 

R + R^+ 

R2 

r + rA 

R2 + I 

R + I 2 

RI + I 

r + rA 

R 2 + r 

Ri + r-^* 

I2 + R 

i+rA 

I 2 + R 
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h formation (reaction rates have units [cm 6 /sec]) 


r + / + rA /a + RI 
I + / + RI / 2 + RI 

r + i + / 2 -^ / 2 + / 2 

i + / + / 2 — —* / 2 + / 2 
/ + / + .ff 2 — / 2 + /i 2 
/ + /* + ii 2 — — + / 2 + fi 2 

Wall reactions (reaction rates have units [cm 3 /sec]) 


/ + I + wall —i+ 7 2 + toail 
f? 2 + / + wall-^-* R + .RJ + wall 

Reaction rate coefficients 

The reaction rate coefficients are assumed to satisfy the following inequalities: 

.476{10“ 16 ) <Qi < 8.4(10 -16 ) 

.730(10“ u ) <Q 2 < 4.94(10“ n ) 

1.23(10“ 18 ) <Q 3 < 11. 1( 10“ 18 ) 

1.57(10“ 16 ) <Q a < 14. 1(10“ 16 ) 

.530(10“ 14 ) <Q 5 < 4.8(10“ 14 ) 


0.9(10“ 13 ) <*=! 
.657(10“ n ) <k 2 
.65(10“ 12 ) <Jb 3 
1.0(10“ 16 ) <k A 
0.33(10“ 11 ) <k 5 
1.0(10“ 17 ) <ke 
1.5(10“ 19 ) <k 7 
0.533(10“ 23 ) <k 8 


< 34.7(10“ 13 ) 

< 8.05(10“ n ) 

< 10.4(10“ 12 ) 

< 9.0(10“ 16 ) 

< 3.0(10“ u ) 

< 10.24(10“ 17 ) 

< 4.5(10“ l9 ) 

< 4.8(10 -23 ) 


1.0(10 -33 ) <c x < 10.2(10“ 33 ) 
1.6(10“ 32 ) <c 2 < 45.0(10“ 32 ) 
4.44(10“ 32 ) <c 3 < 14.4(10“ 32 ) 
2.92(10 -30 ) <c 4 < 4.94(10“ 30 ) 
3.6(10“ 31 ) <c 5 < 6.0(10“ 31 ) 
1.35(10“ 32 ) <c 6 < 2.25(10“ 32 ) 


0.33(10“ 12 ) <vi < 3.0(10“ 12 ) 
0.33(l0“ n ) <V2 < 3.0(10“ u ) 
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Other coefficients 

In the computer program we let: Q y denote the quantum yield; fi, f 2 denote the photodisaociation 
rates which are dependent upon the solar simulator concentration co; <r denotes the stimulated emission 
cross section and p denotes the photon density in the optical cavity. We use the approximations 

fx = 3.04(10" 3 )c o , fa = 3.38(10 _2 )c 0 

The pumping is assumed to occur over an interval 0 < z < where L is the length of the tube. Various 
assumptions can be incorporated concerning the pumping intensity. Currently, the program assumes that 
the maximun pumping intensity occurs at Zqi — 2 z L- 
Differential equations for chemical kinetics 

Using the notation [A] to denote the concentration of species A in units of cm -3 , the differential 
equations defining the chemical kinetics for the iodine lase can be expressed by the following set of coupled 
nonlinear differential equations involving the concentrations [ft/], [.ft], [ft^] , [/ 2 ], [/*], [/]. 

= ki[R]\r\ + * 3 [ft][/] + fc 5 [ft][/ 2 ] - k 7 [I-}[RI\ - k 4 [ft](ft/] 

- * 6 [ft][ft/] - f x [ ft/] + «a[fta][/] - *s[W 

^ = fi m - ki [ft][/*l - *a[ft)[/] - 2Ar 3 [ft] 2 - k 4 [RI\[R\ 

- Jfc 6 [ft/][ft) - * 6 [ft][/ a ] + «*[fta][/] + * 7 [ft/][r] + k s [I\[RI) 

® + kt[Rim + - «*[ w 

® = CifJWlirJI/1 + c 2 [ft/][/] 2 + c 3 [/ 3 ][/*][/] + c 4 [/ 2 ][/] 2 

- fa[/a] + *?[*/][/•] " Wal + «i [I]' 2 

Cs [/] 2 [ft 3 ] + A: 8 [ft/][/] + c 6 [/][rj[ft2] 

® = Q,fi[*/] + fa[/a| - kmn ~ Qa[/a][r ] 

- *°p{\n - 1[/]) + ke[R][Ri] - g 3 [ft]jr] - <? 4 [ft2][r| 

- Qs[r}[I} - fc7[ft/]|/*l - - Cx [ft/] \r\ [/] 

-C 3 [/a] [/•][/]-«! [ft/] [/*] 

^ = &[/a] + <?x[ft/][/*j + Q 2 [/ 2 ][/*] - 2c 5 [/] 3 [ft 3 ] - fc»[/][ft/] 

+ wp(|rj - i[/]) - Cl [ft/][/-][/] - 2c 2 [ft/][/] 2 - c 3 [/ 2 ][/*][/] 

+ 2c 4 (/ 2 ][/] 3 - k 3 [R}[I\ + fc 4 [ft/][ft] + Q 3 (/*][ft] + g 4 [/*][fta] 

+ Q 6 [/*](/] + *s[ftHJa] - va[ft 2 ][/] - 2vi[/] 2 - c 6 [ft 2 ][/*][/] 

In the above differential equations we use the material derivative 

M = M + Mil = M + M 

dt dt dz dt dt dz W 

where ^7 = w is the flow rate in the axial direction. 

The above system of nonlinear differential equations conserves the masses of the species involved in the 
reactions and for steady state operation at any point z we have the immediate integrals 

[ft/] + [ft] + 2[ft 2 ] = constant 
[ft/] + 2[/ 2 ] + [/*] + [/] = constant 
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Photon density 

For the light flux density of lasing photon p we let p = p + + p_ where p + = p+(z, t) denotes the photon 
density propagating in the positive z direction and p~ = p~ ( z , t) denotes the photon density propagating in 
the negative z direction. The differential equations for these photon densities are given by 

l a -§r +d -t-w+u r 

i-iW) 

where c is the speed of light in the optical medium. In the above equations cr\I m ] p+ is the amplification factor 
resulting from population of the upper lasing level of the active medium and — |a[/]p + is the decrease in 
photon density due to population of the lower lasing level. 

Compressible flow 

The effects of fluid density variation as a function of distance z along the tube is considered. Also the 
pressure and temperature of the flow medium are calculated as a function of distance z and incorporated 
into the computer model by including the following equations: 

(i) An equation of state: 

For P the pressure in the gas, T the absolute temperature, and rj the gas density, we assume an equation 
of state for a perfect gas 

P = r\RT 

where R is the universal gas constant. 

(ii) Continuity equation (conservation of mass): 

The continuity equation is expressed 

+ div (ri?) = 0 


where q is the gas density and V is the gas velocity. For steady state conditions and V = uk the flow in the 
axial direction, the continuity equation reduces to 


-(^)=° 


which implies that 

t)u = constant. 

(iii) Momentum equation 

The momentum equation for a control volume having a mass rfdr where dr is an element of volume, is 
given by 

M = JJJ Vridr. 

Using Newton’s second law we have 


F 


DM 

Dt 


D_ 

Dt 


fjf Vr,dr 


where ^ is the material derivative. We have that 




where the surface integral term above represents the efflux of momentum through the control volume and the 
volume integral term represents the change in momentum inside the control volume. Changing the surface 
integral to a volume integral by using the Gauss divergence theorem 


U V(r,V da) = fff[V lVv\ 
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and letting F = fff f dr where / is the average force per unit volume, we obtain 

F - IIJ fir - fff + f^)h] dr 

where W = u> 2 kk is a dyadic and / = — V P is the average force per unit volume which is due to the fluid 
pressure P. This equation implies that 


ap 

dz 


d , \ d . 2 \ 

-(rH + -(^ ). 


Using the result = ci = a constant, the steady state form of the above gives us 


dP du> 
dz Cl dz 


and an integration gives 


P + Cl U/ = C 2 


where c 2 is a constant of integration. 

Energy equation 

In terms of the specific enthalpy h per unit mass, the energy equation for the gas flow can be expressed 



DP 

Dt 


+ kV 2 T + q 


where P is the pressure, T is the absolute temperature, k is the thermal conductivity, and q = q(z) is the 
raditation heat flux. In one dimension, the energy equation can be expressed 


dh 

'dt 


dh 


dP ap <PT 

!t +u> -a; + K d^ + q - 


For C p the specific heat at constant temperature and C v the specific heat at constant volume, we can write 
h = C P T and C p — C„ = R. This gives us the steady state equation 


r ui u J. 

rv{C v (T) + R}—=u—+ K -^+q 

where we have neglected the effects of viscosity. If we also neglect the effects of the thermal conductivity the 
above equation reduces to 

ci[C„(r) + R]dT + ucxdw + qdz = 0 


Using the empirical model 


C V (T) = a„ exp(f) v [T - 300)), 298.15 < T < 500 

where a„ = 147.23 and = 0.0012 sire constants^ the above equation can be integrated to obtain 

Ci— - exp (f3 v (T — 300)) + ciR(T — 300) + -ciw 2 — Q = c 4 


where 

Q = J 9(^) dz 

and C 4 is a constant of integration. 
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Computer program 

The above equations can be found in the computer program listed in the Appendix A. This program 
can be described as follows: 

Main program 

This assigns values to all constants and parameters and then guesses at an initial photon density. The 
equations are then integrated from 0 to £ and the boundary conditions (discussed in an earlier report) are 
checked to see if they are satisfied. If they are not satisfied then an iterative scheme is employed to find the 
initial photon density which satisfies the boundary conditions. When the correct initial photon density is 
used the results of the computations are printed out. 

Subroutine GRAPHS 

Produces graphical output for each of the species concentration as a function of distance z in the axial 
direction. 

Subroutine PFLOW 

This subroutine calculates certain parameters needed in subroutine FLOW. These parameters are stored 
in common BLK10. 

Subroutine FLOW 

This subroutine calculates the temperature T, pressure P, flow rate W, density rj, Pressure in torr, as a 
function of axial distance z. 

Subroutine ARREN 

This subroutine calculates how some of the rate coefficients change with temperature where we have 
assumed various arrenhius expressions for the different rate coefficients. Other rate coefficients are held 
constant. 

Subroutine CHSIl 

This subroutine calculates £1 as a function of z. 

Subroutine CHSI2 

This subroutine calculates ( 2^1 function of z. 

Subroutine COEFFS 

This subroutine calculates the various constants and rate coefficients needed for execution of the pro- 
gram. 

Subroutine VELOC 

This subroutine calculates the velocity w as a function of tube radius. Various assumed flow patterns 
can be assumed. Current version assumes a parabolic flow profile with the velocity of the gas going to zero 
at the tube walls. 

Subroutine FUN 

This subroutine calculates the functions occuring on the right hand side of the differential equations to 
be solved. 

Subroutine SIGMA 

This subroutine calculates the absorption cross section a. 

Subroutine INTEG 

This subroutine integrates the differential equations from 0 to z using a 7th order Runge-Kutta-Fehlberg 
variable step size integration routine. 
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